Invariance of density correlations with charge density in polyelectrolyte solutions 
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We present a theory for the equilibrium structure of polyelectrolyte solutions. The main element 
is a simple, new optimization scheme that allows theories such as the random phase approximation 
(RPA) to handle the harsh repulsive forces present in such systems. Comparison is made with 
data from recent neutron scattering experiments of randomly charged, hydrophilic polymers in 
salt-free, semi-dilute solution at various charge densities. Models with varying degrees of realism 
are examined. The usual explanation of the invariance observed at high charge density has been 
counterion condensation. However, when polymer-polymer correlations are treated properly, we find 
that modeling polymer-counterion correlations at the level of Debye-Hiickel theory is sufficient. 
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Polyelectrolytes are polymers with ionizable groups 
that dissociate in polar solvent leaving ions bound to the 
chain and counterions free in solution. The study of poly- 
electrolytes has traditionally been driven by their om- 
nipresence in living things [l| , but there has been rapidly 
growing use of and interest in synthetic ones for techno- 
logical applications, typically at the nano-scale level 0- 
The theoretical study of them is one of the most chal- 
lenging in the field of complex fluids, however, because 
of the strong repulsions and attractions present in such 
systems 0- 

In this Letter, we present a scheme that allows a poten- 
tially large class of theories to deal properly with these 
strong long-range interactions. We demonstrate the util- 
ity of this approach by comparing with some recent in- 
triguing experiments of Nishida, Kaji and Kanaya (NKK) 
0, and Essafi, Lafuma and Williams (ELW) Counte- 
rion condensation [6j commonly is invoked as the cause 
of the invariance at high charge density observed in these 
experiments. While there is some ambiguity in the liter- 
ature about what precisely is counterion condensation, it 
is generally regarded that the mechanism is not captured 
in any theory in which ion-counterion correlations are 
included at the level of Debye-Hiickel theory 0, Q . How- 
ever, we show in the theory described below that Debye- 
Hiickel modeling of ion-counterion correlations can itself 
produce an invariance at high charge density. 

There are few theories that can handle the harsh re- 
pulsions and attractions in polyelectrolyte solutions ||. 
Further, there is currently no theory that provides re- 
liable predictions for liquid structure in the semi-dilute 
regime — which is the region of most experimental inter- 
est. The theory of Laria, Wu and Chandler has been 
applied successfully to examine density correlations there 
as a function of number density [lOj , but performs poorly 
at large, experimentally relevant, interaction energies and 
charge densities 11]. The theory of Donley, Rajasekaran 
and Liu (DRL) seems to perform properly at large inter- 
action energies and charge densities, but presently only 



for purely repulsive systems [1 2 . Despite its simple form, 
the random phase approximation (RPA) can in some in- 
stances give reasonable predictions for the free energy 
of polyelectrolyte solutions |l3|. Here, we assume that 
theories such as the RPA are also adequate for density 
correlations given a suitable reinterpretion of the inter- 
action potential that appears in them |l4|. Let 

SMkM'k'(r) = ([pMfc(r) - PMk] [pM'fc'( ) ~ PM'k']) (1) 

be the density-density correlation function between a k 
type monomer (or site) on a molecule of type M and a 
k! type monomer on a molecule of type M' a distance 
r = |r| apart. Here, pMfc(r) is the microscopic density of 
fe-type monomers on molecules of type M at position r 
with average value puk- The brackets denote a thermo- 
dynamic average. The Fourier transform of SMkM'k' {r) 
is the partial structure factor and links the theory to scat- 
tering experiments. The correlation function can be sep- 
arated into intra— and inter— molecular pieces, FlMkk' (f) 
and H M kM'k'(r): 



SMkM'k'(r) = &Mkk' (t)SmM' + HukM'k' M , 



(2) 



where Smm' is the Kronecker delta. In this work, the 
molecular structure function Q Mkk' ( r ) is assumed known. 
The intermolecular correlation function i?MfeM'fc'( r ) = 
PMkPM'k'[gMkM'k'(r) - l], where g M kM'k'{r) is the ra- 
dial distribution function. Thus, to compare with exper- 
iment we need a form for this function in terms of the 
potentials and molecular structures. Consider for sim- 
plicity a system where all the interactions can be decom- 
posed as pair potentials UMkM'k'{ r ) between sites. Here, 
we take these potentials to be hard-core for distances r 
less than some range (JMkM'k'i and Coulombic outside. 
Consider also a class of theories for which gukM'k'^f) 
is a function only of the interactions {u} and the total 
molecular structure functions {Q}. Theories in this class 
include the RPA and possible higher loop improvements, 
and the polymer version of the theory of Chandler, 
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Silbey and Ladanyi |lq|. One failure of theories such as 
the RPA is that they predict a negative gMkM'k'{r) for 
short distances if the interaction u MkM > k >(r) is strongly 
repulsive even though the radial distribution function is 
intrinsically non-negative. This shortcoming is thought 
to be one reason that RPA underestimates the strength 
of the dependence of 

Qmax with p m in the semi-dilute 

regime. 

For short-range repulsive, e.g., hard-core, potentials 
one remedy to this deficiency was proposed by Ander- 
sen and Chandler(AC) 14]. Their idea was to replace 
the true potential u by a pseudo or optimized one, 
UMkM'k'(r)- For distances r < a M kM'k', u M kM'k'(r) 
takes a value such that the hard-core condition is satis- 
fied, i.e., guikM'k'ir) = 0. Outside the core, u M kM'k'{r) 
can be set equal to the original potential u which is long- 
ranged and presumably weak. In this way, higher order 
diagrams neglected in the theory are summed in an ap- 
proximate manner to enforce the core condition. When 
applied to the RPA, this optimization scheme is a gener- 
alization of the MSA closure of the Orstein-Zernike equa- 
tion 0] and a closely related scheme has been shown to 
be diagramatically proper for CSL theory [l^. This op- 
timized RPA or ORPA, usually with long-range poten- 
tials absent, is more popularly known as RISM for small 
molecules and PRISM for polymers 0. 

Unfortunately this AC optimization is not very useful 
for polyelectrolytes. One reason is that the contact en- 
ergy between ion and countcrion seems to be less impor- 
tant than the interaction energy of the counterion with 
the whole chain. Hence, enforcing the hard-core behav- 
ior between opposite charges matters less than modeling 
correctly the chain structure, e.g., determining the bond 
length b. A more important reason is that the repulsion 
between like-charged polymers is usually very large. This 
causes the radial distribution function to be effectively 
zero out to a distance that scales with the Debye-Hiickel 
screening length for semi-dilute densities which is usu- 
ally much larger than the hard-core distance 0, 0, • 
Hence, enforcing the core condition produces little im- 
provement in the theory for long polyelectrolytes. 

An alternative scheme is to optimize the range of the 
pseudo-hard-core portion of u, and not just its ampli- 
tude. By this we mean that if g(r) is nearly zero out 
to some distance, then it makes no difference whether 
this exclusion zone is caused by a hard-core potential 
or a Coulombic one. This effective hard-core diameter 
a MkM>k> 1S determined by requiring that gMkM'k'(r) be 
non-negative everywhere (not just zero inside the core). 
The closure to the theory then is: 

gMkM'k'(r)= 0, r << J M f kM>k>> (o\ 

UMkM'k'{r) = UMkM>k>{r), r > <T^ IkM , k ,, 

and the range cr^{ M , fc , of the hard-core interaction is cho- 
sen to have the smallest value such that gMkM'k'if) > 



for r > a MkM'k' subject to the constraint that 

a MkM'k' ^ AMkM'k'- With guesses for the &MkM'k' , ^'> 
one can solve the RPA (or similar theory) numerically in 
a manner very similar to what is done in PRISM. One 
then adjusts the c^ffeM'*;' by iteration. Eq. © is the 
primary result of this paper. 

For simplicity, we show results for the range opti- 
mized version of RPA, denoted as RO-RPA. We com- 
pare with scattering data of NKK and ELW. In these 
experiments, the molecules were hydrophilic, sulfonated 
vinyl polymers that were linear and randomly charged, 
with monovalent ions and counterfoils. Small angle X-ray 
(SAXS) and neutron (SANS) scattering methods were 
used to examine the equlibrium structure of polyelec- 
trolytes and counterions in their salt-free solution, water. 
Previous experiments on polyelectrolytes have examined 
how density correlations in the liquid change with poly- 
mer monomer number density p m . Here, though, p m was 
held constant at a semi-dilute or near semi-dilute density 
and the per chain average fraction of charged monomers 
/ was varied. Batches of polymers with different values 
of / were created and the scattered intensity I(q) as a 
function of wavevector q was measured for each. NKK 
found that at small /, the peak position of I(q), q m axi 
obeyed a power law, but at larger / ~ 0.4 appeared 
to reach an asymptote. ELW extended these measure- 
ments to higher / and found that q ma x was effectively 
constant for / > 0.4. SANS measurements also allowed 
them to extract the polymer monomer-monomer struc- 
ture factor S mm (q). They found that at least for the 
range 0.55 < / < 0.81, S mm (q) was invariant for the 
wavevectors measured. A picture that ELW offer is that 
as one increases / above 0.4, the system acts as if the 
charges on the chain are renormalized such that the ef- 
fective / is constant. Counterion condensation is usually 
given as the explanation. 

The neutron scattering data of ELW on poly- 
AMAMPS is sufficiently well characterized to allow quan- 
titative comparison via S mm {q) 0. We show results for 
Smm{q) for three models, the "minimal", "primitive" and 
"two-state". In all models the polymer chains are lin- 
ear and rod-like. The chains are uniformly charged and 
all monomers are identical so that the monomer valency 
Z m = fi which can be varied continuously from to 1. 
The strength of the system interactions is characterized 
partly by the ratio of the Bjerrum length lg — e 2 / (eksT) 
to the polymer bond length b. Here, e, e and fc^T are 
the electron charge, solvent (water) dielectric constant 
and thermal energy, respectively. In the models here, 
the effect of the salt-free solvent enters only through e. 
Both monomer and counterion are spherical and we set 
their diameters, a, equal to each other and equal to the 
bond length b. For ELW, l B ~ 7.lA, b « 2.5 A and 
p m <j 3 « 3 x 10~ 3 . The counterion density p c was de- 
termined by charge neutrality, Z m p m + Z c p c = and 
counterion valency Z c = — 1. The chain length N was set 




FIG. 1: Scaled peak wavevector q m ax<? of the monomer- 
monomer partial structure factor S m m{q) as a function of 
the average fraction of charged monomers per chain / for 
?b/6 = 2.85 and p m cr 3 = 3 x 10~ 3 . The meaning of the curves 
and symbols is shown in the figure legend. 



FIG. 2: Scaled peak height of the monomer-monomer partial 
structure factor Smm{qmax) / pm as a function of the average 
fraction of charged monomers per chain /. The conditions 
are the same as in Figure 1. The meaning of the curves and 
symbols is shown in the figure legend. 



to 500. In the minimal model the system consists only 
of linear polymers. Monomers on the same chain or dif- 
ferent chains interact with a potential which is hard-core 
for r < a. Outside the core the potential has a screened, 
Debye-Hiickel form: u mm (r) = Zm ^ B e~ Kr , where the in- 
verse screening length k = y/ ' AttZ^IbPc- In the primi- 
tive model, the system consists of linear polymers and 
counterions. Polymer monomers and counterions inter- 
act with a potential which is also hard-core for r < a. 

Outside the core the potential has a bare Coulomb form: 

z ■ z i 

u ij ( r ) = ' r 3 B , where i or j denote monomers or coun- 
terions. 

Comparison of the predictions of RO-RPA in the prim- 
itive model with simulation data [2l| shows that RO- 
RPA tends to underestimate the strength of polymer- 
counterion correlations 01 The two-state model is a 
common way to improve these correlations 0, llflj . In 
this model counterions are divided into two species, free 
and condensed. In one version of the model, free ones 
are counterions as in the primitive model, and condensed 
ones are bound to the surface of the chain, but able to 
translate along its length. Each chain has the same num- 
ber of condensed counterions, N c . An average N c is 
determined by constructing a free energy for the sys- 
tem and then minimizing it with respect to N c . For 
a given N c we computed the intramolecular structure 
functions fijwfcfe' ( r ) f° r the polymer-condensed counte- 
rion "molecule" by a single chain Monte Carlo simulation. 
Here, the intramolecular effective interaction between 
chain monomers and condensed counterions was taken to 
have a screened Debye-Hiickel form similar to above, but 



the inverse screening length n = ^/AirlBiZ^pm + Z^p c ) 
was due to both ions and counterions. We calculated 
the free energy using the "charging" formula [l^, as- 
suming that the molecular structure was constant during 
the charging. This constancy is not correct obviously; 
however, the goal here is not a value for the free energy, 
but the position of its minimum with respect to N c . We 
expect for the minimum that this approximation is a rea- 
sonable one. 

Figures 1 and 2 show results for the peak position q max 
and peak height S mm (q max ) of the monomer- monomer 
partial structure factor. The solid, dashed and dashed- 
dotted lines correspond to results of RO-RPA using the 
minimal, primitive and two-state models, respectively. 
Results from the RPA are also shown as dotted lines. 
Diamonds are data of ELW. At small /, both quantities 
seem to obey power laws: q max ~ f a and S mm (q max ) ~ 
/-'>. All RO-RPA models predict a « 2/3 and ry « 0.9. 
For flexible hydrophilic polyelectrolytes a is expected to 
be 1/3 0, I23. The difference is due probably to our 
modeling the chains as rods while they are predicted to 
be rod-like only for a restricted range of lengthscales (2^ • 
We know of no data for 77 for polyelectrolytes. At large 
/ > 0.4, it is clear that range optimization completely 
changes the character of the RPA theory. In contrast 
to the RPA, RO-RPA in the minimal model predicts 
that both quantities become pretty much constant, in 
agreement with NKK and ELW @|. Note that adding 
explicit counterions and condensed counterions produces 
only moderate improvement in the theoretical trends |23| . 
As a check on the validity of the RO-RPA theory, we 
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also show results (dot-dot-dashed lines) of the DRL the- 
ory |l2| in the minimal model. DRL is an approximation 
to the "two-chain" equation for g(r) [9|,|24j. As can be 
seen, this theory also produces an invariance at large / 
and agrees almost quantitatively with RO-RPA at all val- 
ues of /. Agreement between RO-RPA and experiment 
for S mm (q max ) is very good with the former predicting 
values less than 5% lower than the latter. Agreement for 
<lmax is less satisfactory with RO-RPA predicting values 
about 35% higher than experiment. The most probable 
cause of this discrepancy is our neglect of chain flexibil- 
ity since the Debye-Hiickel screening length £ ~ 5<r at 
large / implying a flexible chain on scales larger than 
that. If the chains were partially collapsed, then the re- 
pulsion and average distance between monomers on dif- 
ferent chains would be larger. Thus q max would decrease 
and S mm (q max ) would increase. On the other hand, if 
the solvent — which is "good" for poly-AMAMPS — were 
included explicitly then S mm (q max ) would decrease. An 
explanation for this invariance with / seen by the theory 
is due to screening: as / increases, the repulsion between 
polymer chains increases; on the other hand, charge neu- 
trality forces more and more counterions into the solu- 
tion which increases the screening between polymers. At 
small / polymer and counterion correlations are weak, 
but at large / they increase such that one reaches a bal- 
ance between polymer-polymer repulsion and counterion 
screening. 

In conclusion, we have presented a new optimization 
scheme that appears to improve substantially the pre- 
dictive power of the RPA theory for the structure of 
polyelectrolyte solutions. We showed that this range- 
optimized RPA in the minimal model yields all trends 
exhibited in the experimental data of ELW at large / 
and is in moderate agreement with that data for the peak 
position and height. The crudeness of our chain model, 
however, precludes us from drawing any firm conclusions 
about the origins of the effects seen in ELW. On the other 
hand, since the minimal model does predict an invariance 
and the comparison is reasonable enough, it is apparent 
that a number of questions need to be answered before 
the cause(s) of this phenomenon is determined. But it is 
interesting that a model of screened Coulomb chains can 
produce effects that would normally be associated with 
more complex mechanisms such as counterion condensa- 
tion. 
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